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X-ray absorption near-edge spectroscopy (XANES) is used to study the N environment in bulk 
GaN and in GaN y Asi_ H epilayers on GaAs (001), for y ~5%. Density-functional optimized struc- 
tures were used to predict XANES via multiple-scattering theory. We obtain striking agreement for 
pure GaN. An alloy model with nitrogen pairs on Ga accurately predicts the threshold energy, the 
width of the XANES "white line", and features above threshold, for the given X-ray polarization. 
The presence of N-pairs may point to a role for molecular N2 in epitaxial growth kinetics. 

PACS numbers: 61.10.Ht,81.15.Hi,71.15.Mb 



Nitrogen- alloyed GaAs, i.e., GaN y Asi_y, is intensely 
studied owing to its technological importance. Even 
very small amounts of N lead to a dramatic y-dependent 
decrease in the bandgap l]. In the ultradilute regime 
(y < 0.01%), localized, single-impurity levels produce 
sharp resonances above the GaAs conduction band min- 
imum. With increasing y, small nitrogen aggregates 
(pairs and higher clusters) give rise to localized states in 
the bandgap, producing sharp photoluminescence (PL) 
lines 0,0- At higher concentrations {y > 0.1 — 0.25%), 
the sharp PL lines merge to a broader line. Anomalous 
bandgap bowing and other effects are observed 0> Q in 
this regime. At these concentrations, a phenomenologi- 
cal band anti-crossing model works well [5J, |6J , since here 
the details of N-incorporation seem to be unimportant. 
Kent and Zunger[4| presented a theoretical study using 
pseudopotentials and large supercells containing various 
N-atom configurations. In reality, a plethora of different 
local atomic environments is possible, since Ga can have 
to 4 first neighbour N atoms; N pairs can have dif- 
ferent orientations or separations; and triplets or higher 
order clusters are possible 0, Q • The change in PL on 
annealing has even been attributed to the removal of N 
interstitials or defect complexes, such as N-As split in- 
terstitials 0,0. 

Several groups have used XANES to study local atomic 
structures in GaNAs epilayers. Lordi et al. 10] mod- 
eled their N XANES data using local densities of states 
(LDOS) obtained from the Vienna Ab Initio Simula- 
tion Package(VASP) [ll|. For GaNAs the partial LDOS 
did not reproduce the experimental "white line" and 
predicted several unobserved oscillatory features above 
the absorption edge. Simpler tight-binding (TB) meth- 
ods have also been examined by other authors [L| . In 
general TB is unreliable except for band-structure type 
calculations^) to which the TB parameters have been 
fitted. 

In this Letter we present a study of N incorporation 



in GaNj,Asi_j, epilayers. Our objectives are: (i) to com- 
pare the local nitrogen environments in the dilute nitride 
and in bulk GaN; (ii) to predict the XANES spectra for 
several local atomic arrangements of N in GaN^Asi-^ 
and compare them with experiment. These calculations 
used the ab-initio code FEFF8 which combines si- 
multaneous real-space full multiple scattering (FMS) cal- 
culations of near-edge X-ray absorption spectra within 
a self-consistent field (SCF) treatment of the electronic 
structure. The microscopic structure of the alloy epilay- 
ers for N or (N,N) pair-configurations were separately 
determined via density-functional theory (DFT) using 
the VASP code 11] within an ultrasoft pseudopotential 
scheme, and the generalized gradient approximation for 
the exchange-correlation potential. The DFT optimized 
structure was used as the input to the XANES calcula- 
tions. The predicted XANES spectra are in very close 
agreement with the experimental results, and support 
the presence of (N,N) pair clusters in this alloy regime 
(y ~ 3 - 6%). 

The single, uncapped GaNj,Asi_j, epilayer used 
here was grown by solid-source MBE as described 
elsewhere |l6|. The actual composition, y = 0.0496, and 
thickness, 333 A, of the sample were obtained from high- 
resolution X-ray diffraction data via dynamical diffrac- 
tion analysis. A thick, metalorganic vapor phase epitaxy- 
grown GaN layer on sapphire (0001) was used as a ref- 
erence. Experimental XANES measurements were per- 
formed at the Canadian Spherical Grating Monochroma- 
tor beamline of the Synchrotron Radiation Center at the 
University of Wisconsin, Madison. The storage ring en- 
ergy was 1 GeV, and a 600 lines/mm grating was used. 
The measurements on GaN and GaNAs epilayers used 
entrance-slit widths of 100 /im and 250 /jm, providing 
energy resolution of ~ 0.24 eV and 0.54 eV, respectively. 
Measurements at the nitrogen K-edge were performed 
using both total electron yield (TEY) and total fluores- 
cence yield (TFY) modes of detection. TEY provides 
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more surface-sensitive information. For GaN, the TEY 
and TFY data were nearly identical, but the TEY data 
exhibited a more intense main peak. For GaNAs, only 
the TFY data are presented here. The samples were out- 
gassed by in- situ heating to 280 °C for ~4 hours in the 
ultra-high vacuum analysis chamber. The effectiveness of 
this surface preparation was evident from the elimination 
of the Oxygen K-edge peak in the absorption spectra. 

Five GaN y Asi_ y models were considered in this study. 
For each structure, state-of-the-art DFT calculations 
within the generalized gradient approximation (GGA) 
were carried out with a 64 atom 2x2x2 supercell of 
zincblende unit cells, with periodic boundary conditions. 
Tetragonal distortions were included by holding the in- 
plane epilayer lattice constant at the GaAs (substrate) 
value. The lattice constant along the growth direction 
and atomic positions were relaxed until the Hellman- 
Feynman forces became negligible, resulting in the lat- 
tice constant along the growth direction being smaller 
than the in-plane value. Each FEFF model was formed 
by "stacking" 27 relaxed DFT supercells to form a 1792- 
atom, cubic supercell with the absorbing-N atom at the 
center. Spherical clusters using only the 295 central 
atoms were used for the FEFF calculations. A single 
N atom in the DFT cell yields the alloy concentration 
y = 1/32. The artificial periodicity of the DFT supercell 
resulted in second N atoms at ~ 11.2 A from the central, 
"quasi-isolated" N absorber. All our models contain this 
artificial periodicity. Using the notation of Ref. 0, we 
also considered (N,N) pairs in 771th nearest neighbour po- 
sitions with rn=l, 2 and 6, corresponding to FCC lattice 
positions (1/2 1/2), (1 0) and (I 1 I). For the m=2 
case, we studied (N,N) pairs parallel (m=2a) or perpen- 
dicular (m=2b, (0 1)) to the substrate, with the X-ray 
absorbing atom at (0 0). Thus, we have one "single N" 
model, and four (N,N) pair models labeled as NNto, i.e., 



Single N, NN1, NN2a, NN2b and NN6 



(1) 



For nitrogen K-edge XANES, large FEFF clusters are 
needed to capture multiple scattering resulting from the 
high density and the long mean free path resulting from 
the very long core-hole lifetime (inverse lifetime = 0.116 
eV). In this work, 27 and 295-atom SCF and FMS clus- 
ters were used and parallel-processed on 32 nodes. 

Fig- HI shows the experimental data for the thick GaN 
reference epilayer obtained at incidence angles of 0° and 
45°, relative to the surface normal. The data show a rich 
fine structure consistent with previous reports 0, 0] . 
These two data sets illustrate the importance of polar- 
ization effects, (cf. Ref. 0). Most notably, the peaks 
near 402.5 eV and 406.8 eV are significantly enhanced at 
45° compared to normal incidence. The former is close to 
the "magic" angle, 54.7° (the Lambert angle), where the 
polarization effects would be completely suppressed . 

The figure also shows the FEFF8 results for the GaN 
N K-edge XANES based on the hexagonal Wurtzite 
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FIG. 1: Nitrogen K-edge XANES for wurtzite GaN compared 
with theory (unpolarized calculation). Incidence angles of 0° 
and 45° were used in single scans at a rate of 10 sec/point, and 
step size of 0.25 eV. The inset shows the excellent agreement 
for incidence at 45°, close to the Lambert angle. 



structure with the unit cell parameters a = 3.186 A, c 
= 5.178 A. Clearly, the theory is very successful in pre- 
dicting the main XANES features, especially at 45°. In 
particular, the edge energy is accurately predicted with- 
out any "fitting" to experiment. The experimental ab- 
sorption edge energy, E — 401.1 ± 0.1 eV was deter- 
mined from the zero of the second derivative of the ab- 
sorption, while the value Eq — 400.93 eV was predicted 
by FEFF. The experimental uncertainty is about half 
the energy step (0.25eV). On normalizing to match the 
height of the main peak, ~405 eV, the theory accurately 
predicts the magnitudes and positions of the peaks near 
402.5, 406.8 and 411 eV, as well as the local minimum 
near 414 eV, the local maximum near 423 eV, and the 
general decay of the absorption with increasing photon 
energy. A determination of the experimental polariza- 
tion was difficult given the irregular shape of the sample 
after cleaving the sapphire substrate. Hence the appro- 
priate comparison is between our unpolarized model and 
the experimental data at 45°. Unlike the theoretical re- 
sults of Ref. flTl our results show no overestimate of the 
peak positions at higher energies. This suggests that the 
multiple-scattering methods of FEFF8 can be used for 
GaNyAsi-y with confidence, if a realistic alloy structure 
is available. 

Fig. |3 shows the experimental XANES for the 
GaN y Asi_j, epilayer. The local environment of the Ni- 
trogen in the alloy is clearly different from the pure GaN. 
The absorption threshold in this case has a single, nar- 
row "white line" , as opposed to the multiple peaks in 
GaN. A smaller, distinct peak is observed at 402.3 eV, 
coincident with the lowest energy peak in pure GaN. The 
top panel shows that the NN1 model accurately predicts 
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FIG. 2: Nitrogen K-edge XANES for a coherently-strained 
GaNAs epilayer (thick line) compared with theoretical models 
listed in Eq.0 The data are single scans, with 135 sec/point, 
and a step size of 0.25 eV. Top: for photons polarized along 
(110). Bottom: unpolarized photons. The XANES for NN2b 
is very similar to NN2a, and is omitted. The theoretical ab- 
sorption asymptote for higher eV is lower than in experiment 
as self-absorption effects are not included in the theory. 



the first and second peaks when the photon polarization 
is included in the theory. The lower panel shows that 
the unpolarized calculation may have some relevance to 
the threshold edge, and for energies above ~ 405 eV. The 
threshold energy was found to be Eq — 399.6±0.1eV. Our 
experimental data are consistent with those in Refs. Ilfj 
andll3l The data of Ref. are similar, although their 
second peak intensity is almost at the "white line" . 

Several properties determined from the DFT calcula- 
tions are summarized in Table [J The DFT predicted 
tetragonal distortions are in agreement with Vegard's 
law, confirming its applicability in this alloy regime. 

For the single-N atom model, the Ga-N bond lengths 
are nearly identical, R = 2.0567A, while the NN1 model 
predicts a distribution of bond lengths, R = 2.072 ± 
0.042A. This is an indication of lattice distortions pro- 
duced by this configuration, and we find that this (N,N) 
pair configuration is not as energy optimal as some 
other models. However, Molecular Beam Epitaxy is a 
non-equilibrium technique hence the MBE-grown epi- 
layer need not have the most energy-optimal N-cluster. 
The objective of the XANES is in fact to find the most 



TABLE I: Growth direction (001) lattice constants a±, and 
average Ga-N bond lengths b with standard deviation <jf,, and 
nearest nitrogen distance R obtained by total energy mini- 
mization for the alloy models (see Eq. 0in the text) studied 
here. N.B., the Ga-N and Ga-As bond lengths in cubic GaN 
and GaAs are 1.9500 and 2.4480 A, respectively. The in-plane 
lattice constant is constrained to the GaAs value, 5.6534 A. 



Model 


a ± (A) 


6(A) 


at (A) 


R (A) 


Single N 


5.5865 


2.0567 


0.0002 


11.1730 


NN1 


5.5131 


2.0722 


0.0416 


3.6908 


NN2a 


5.5108 


2.0504 


0.0010 


5.6530 


NN2b 


5.5195 


2.0510 


0.0007 


5.5191 


NN6 


5.5087 


2.0468 


0.0010 


9.7066 



prevalent N cluster in the alloy sample prepared via 
MBE. 

Fig. [21 bottom panel illustrates the spectral changes 
due to the local N environment in the singlet and pair 
models, not including polarization effects. In this figure 
the FEFF8 spectra were shifted (by an amount smaller 
than the experimental stepsize 0.25eV) to match the ex- 
perimental edge energy. The threshold energy E t h is 
almost identical for the different models, and hence in- 
dependent of the local bonding configurations and the 
number of N atoms in these models. The shift in the 
E t h between GaN and GaN y Asi_ y (-1.45 eV) is thus a 
general feature of N alloying in GaAs, and is related to 
long-range effects and the larger lattice constant of the 
GaNyAsi-y epilayer compared with GaN. This is con- 
sistent with the concept of nitrogen impurities forming 
"perturbed host states" in GaAs [j]. 

The apparent discrepancy between the calculated in- 
tensities at higher energies compared with the experi- 
mental results is due to a reduced white line intensity in 
our fluorescence data due to self-absorption effects; our 
normalization using this lower peak intensity causes the 
apparent increase in the absorption at higher energies. 
Otherwise, the "white line" is reproduced extremely well 
by FEFF, and is not affected by the local nitrogen envi- 
ronment. The XANES white line is often understood as 
arising from a high density of final states or exciton effects 
due to the Coulomb interaction between the photoelec- 
tron and the core hole . Nodwell et al. suggested that 
the white line is due to a N-related resonance in the con- 
duction band and not to an excitonic bound state [l3j |. 
However, core hole effects were included in our FEFF 
calculations and provide an excellent description of the 
data. 

The NN6 model and the single-N atom model have 
very similar spectra (Fig. bottom panel), indicating 
that the phase interference between the two atoms with 
the large NN6 separation is sufficiently small as to have 
no effect on the XANES, in agreement with the negligible 
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interaction energy for m = 6 pairs (see also Ref. |4J). In 
the models where the second N atom is moved closer to 
the N absorber, the interaction as well as interference 
of scattered amplitudes becomes more significant. For 
the NN2 models interactions produce a slight increase in 
the intensity of the high-energy shoulder near the white 
line, but the rest of the NN2a and NN2b spectra are 
similar to the single and NN6 models. A more dramatic 
effect is observed for the NN1 model. The intensity of the 
shoulder is increased and the shoulder is shifted closer to 
the small peak at 402.3 eV of the experimental data. 
All models predict a minimum near 413 eV, well-aligned 
with the falling edge of the final experimental peak. The 
predicted features between 413 eV and about 419 eV are 
common for all models and similar to those in GaN in this 
energy range. These spectral features are due mostly to 
Ga-N bonds, traditionally studied by EXAFS. 

We return to the polarized photon calculation of Fig.|2] 
The absorption cross-section depends on the components 
of atomic bonds along the direction of the electric field 
polarization. For the NN1 pair model, the NN1 neigh- 
bour is located at (2.6438A, 0.000A, 2.5753A), so the 
(110) and (110) components are approximately equal. 
The agreement between theory and experiment is strong 
evidence that NN1 pairs are responsible for the peak at 
402.3 eV in the experimental data. We also note that 
the 402.3 eV peak is not shown by the other models and 
is not well-resolved in several triplet models (see Ref. 0) 
that we examined, but have not reported in detail in 
this paper. If the tetragonal distortion (coherent epi- 
layers) arising from fixing the in-plane lattice constant 
to the subtrate is relaxed, the resulting "annealed" NN1 
model loses the 401.4 eV minimum (fig. 2 top panel shows 
only the coherent epilayer results) and the 402.3 eV peak 
moves to lower energies and becomes a shoulder. This 
type of peak narrowing is expected since the annealing 
reduces the bond-length distribution. In fact, the width 
of the 402.3 eV peak in the experimental data (fig. 2, top) 
as compared to the theory suggests that the samples have 
more disorder than in the optimal tetragonal NN1 model 
used in the theory. 

Nodwell et al. noted that their XANES data were 
the same for samples with different N compositions |l3j . 
This seems reasonable since the observed fine structure 
contains major contributions from microscopic N aggre- 
gates, such as nitrogen first-neighbour FCC pairs, and 
these structures are formed even at very low N content. 
Ref.0 considered several triplet structures, includin g th e 
Cs symmetry triplet proposed by Gil and Mariette |2Cj . 
which was found to be consistent with GaNP data. The 
cluster formation does not result solely from the natural, 
statistical distribution of N atoms. The low miscibility 
of N in GaAs has been suggested as a cause of locally 
ordered structures [2lt . An interesting possibility for 
the origin of close N,N pairs is that they are favoured by 
the kinetics of the epitaxial growth process. In fact, it 



has recently been proposed that molecular nitro gen in an 
excited state is a dominant agent in GaN MBE.|23j 

In conclusion, high-resolution XANES of a coherently- 
strained GaN y Asi_ a epilayer on GaAs(OOl) and a bulk 
GaN film were presented. For pure GaN, our theoret- 
ical spectrum agrees extremely well with the measured 
spectrum and demonstrates the effectiveness of the full 
multiple scattering approach within the self-consistent 
field calculations of the Fermi energy, and the impor- 
tance of the synchrotron polarization. For the strained 
GaN y Asi_ y epilayer on GaAs(OOl), the calculated spec- 
tra exhibited a strong sensitivity to the local N environ- 
ment and we show that nitrogen FCC pairs attached to 
Ga captures the observed XANES. 
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